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Summary:  General  variational  theorems,  for  the  rate  problem  of  classical 
elasto-plasticity  at  finite  strains,  in  both  Updated  Lagrangean  (UL)  and 

Total  Lagrangean  (TL)  rate  forms,  and  in  terms  of  alternate  measures  of  stress- 
rate  and  conjugate  strain-rates,  are  critically  studied  from  the  point  of 
view  of  their  application.  Attention  is  primarily  focussed  on  the  derivation 
of 'Consistent  complementary  energy  rate  principles,  which  could* form  the  • 
basis  of  consistent  and  rational  assumed  stress  type  finite  element  methods; 
and  two  such  principles,  in  both  to  UL  and  TL  forms,  are  newly  stated.  Systema- 


tic procedures  to  exploit  these  new  principles  in  the  context  of  a finite 
element  method  are  also  discussed.  Also  certain  general  modified  variational 
theorems,  to  enable  an  accurate  numerical  treatment  of  near  incompressible 
behaviour  at  large  plastic  strains,  are  discussed. 

I.  Introduction 

The  advent  of  high  speed  digital  computers  and  powerful  numerical  methods, 
such  as  the  finite  element  methods,  in  the  past  two  decades  or  so,  have  great- 
ly expanded  the  scope  of  application  of  nonlinear  theories  of  solid  continua 
to  practical  problems  in  engineering.  In  the  formulation  of  such  numerical 
methods  as  the  finite  element  methods,  for  problems  of  nonlinear  solid 
mechanics,  variational  theorems  (and  their  generalizations  to  account  for 

J 


discontinuities  at  interelement  boundaries  in  a finite  element  assembly) 

have  played  a central  role.  [See,  for  instance,  WASHIZU  (1975),  NEMAT -NASSER 

(1974),  ATLURI  (1975),  and  ATLURI  and  MURAKAWA  (1977),  for  a discussion  of 

finite  element  formulations  in  nonlinear  elasticity]. 

Rigorous  and  consistent  formulations  for  numerical  analysis  of  large 

strain  elasto-plastic  problems  have  become  necessary  due  to  the  increased 

interest,  in  recent  years,  in  analyzing  problems  such  as  metal-forming 

processes,  ductile  fracture  initiation  and  stable  crack  growth 

in  cracked  bodies,  etc.  Indeed  several  such  formulations,  and  applications 

of  the  same,  have  appeared  in  recent  literature.  Among  these  can  be  cited 

the  works  of:  HIBBIT,  MARCAL,  and  RICE  (1970),  who  use  a total  Lagrangean 

(TL)  rate  formulation  [wherein  a fixed  reference  frame  is  used] ; NEEDLEMAN 

(1972),  NEEDLEMAN  and  TVERGAARD  (1977),  and  HUTCHINSON  (1973)  who  also  use  a 

...»••  * • * 

.a  ■%  9 # # 

TL  rate  formulation,  but  with  convected  coordinates;  YAMADA,  HIRAKAWA,  and 
WIFI  (1977)  who  use  an  Updated  Lagrangean  (UL)  rate  formulation  [wherein 
the  current  configuration  is  used  as  a reference  for  the  subsequent  step]; 
OSIAS  (1972)  who  also  uses  an  UL  scheme, which , due  to  the  use  of  an  elastic- 
plastic  rate  constitutive  law  that  does  not  admit  to  a potential,  leads  to 
non-symmetric  stiffness  matrices  through  a Galerkin  Scheme;  McMEEKING  and 
RICE  (1975)  who  also  use  a UL  scheme  which,  through  the  use  of  a rate  consti- 
tutive law  with  a potential,  leads  to  symmetric  stiffnesses;  and  NEMAT-NASSER 
and  TAYA  (1976)  whose  formulation  represents  a modification  to  that  of  Mc- 
MEEKING and  RICE  (1975)  to  improve  the  accuracy  in  the  case  of  large  deforma- 
tion of  compressible  materials.  All  the  above  cited  works  employ  a classical 
rate- independent  elasto-plastic  theory,  as  generalized  by  HILL  (1959).  It 
should  also  be  noted  that  all  the  above  finite  element  rate  formulations 
are  based  on  the  principle  of  virtual  work  in  rate  form,  as  first  stated  by 
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HILL  (1959);  and  thus  all  the  above  finite  element  schemes  are  based  on 
assumed  displacements  that  are  compatible  at  interelement  boundaries. 

However,  to  the  best  of  the  author's  knowledge,  no  studies  concerning 
the  convergence  of  the  assumed  displacement  finite  element  methods  of  the 
above  cited  type,  for  the  rate  problems  of  classical  elasto-plasticity , exist 
in  the  literature.  Even  in  the  somewhat  simpler  problem  of  finite  elasticity, 
studies  of  convergence  of  finite  element  methods,  based  on  potential  energy 
principles,  are  just  beginning  to  emerge  [ODEN  (1978)].  From  this  stand- 
point, as  well  as  from  that  of  possibly  studying  solution  bounds,  it  is  of 
interest  to  consider  consistent  formulations  of  numerical  (finite  element) 
methods  based  on  complementary  energy  principles  for  the  rate  problem  of 
finite  strain  elasto-plasticity. 

Another  important  question  in  numerical  schemes  for  elastic-plastic  flow 

• • • * • • ' •* 

at  large  strains  is  how  to  deal  with  the  effectively  incompressible  behaviour 
at  such  magnitudes  of  strain.  It  is  well-known  that  numerical  schemes  based 
directly  on  the  principle  of  virtual  work  fail  in  the  limit  of  incompressibility 
unless  the  mean  stress  is  introduced  as  an  additional  variable  in  the  formu- 
lation. Such  formulations,  which  are  essentially  variations  of  the  well- 
known  HELLINGER  (1914)  - REISSNER  (1950)  theorem,  were  introduced  for  nearly 
or  precisely  incompressible  linear  elastic  materials  by  HERRMANN  (1965), 

TAYLOR,  PISTER,  and  HERRMANN,  (1968),  and  by  KEY  (1969).  To  improve  the 
numerical  accuracy  in  the  near- incompressible  cas*e,  NAGTEGAAL,  PARKS,  and 
RICE  (1974,  Appendix  2)  modify  their  UL  rate  potential  energy  formulation  for 
elasto-plasticity,  in  a way  analogous  to  that  of  KEY  (1969),  except,  instead 
of  the  mean  pressure  as  in  KEY  (1969),  they  use  the  dilatational  strain  rate 
as  an  independent  variable.  As  a consequence,  even  though  the  formulation  of 
KEY  is  valid  for  both  nearly  and  precisely  incompressible  cases,  the  formula- 
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tion  of  NAGTEGAAL  et.al  ceases  to  be  valid  in  the  case  of  precise  incompressi- 
bility. However,  due  to  the  inherent  nature  of  the  complementary  energy 
principle  (with  assumed  stresses  as  variables),  it  is  much  easier  to  treat 
situations  of  near  or  precise  incompress iblity  when  finite  element  schemes 
based  on  a complementary  energy  principle  are  used.  The  works  of  TONG  (1969), 
and  PIAN  and  LEE  (1976)  in  solving  problems  of  linear  elastic  near- incompres- 
sible solids;  and  that  of  MURAKAWA  and  ATLURI  (1978b)  in  solving  finite 
strain  problems  of  incompressible  nonlinear  elastic  solids,  all  of  which  are 
based  on  appropriate  complementary  energy  principles,  tend  to  support  the 
above  view.  Moreover,  as  is  well-known,  better  solutions  for  stresses  are 
obtained  from  numerical  schemes  based  on  complementary  energy  principles  than 
from  those  based  on  potential  energy  principles  (wherein,  the  stress  solution 
is  obtained  by  differentiation  of  the  numerical  solution  for  displacements, 
which  results  in  a loss  of  accuracy  for  stresses). 

Also,  as  noted  by  MASUR  and  POPELAR  (1976),  the  complementary  energy 
approach  holds  a considerable  promise  for  applications  to  buckling  problems, 
wherein,  an  approximation  to  the  stress  state  before  buckling  is  often 
possible  even  when  displacements  remain  unknown. 

For  the  above  reasons,  special  emphasis  is  placed  in  the  present  paper 
on  the  study  of  the  existence  of  consistent  complementary  energy  rate  theorems 
for  the  rate  problem  of  finite  strain  classical,  rate- independent,  elasto- 
plasticity.  Both  the  types  of  formulations,  viz.,  the  Total  Lagrangean  as 
well  as  the  Updated  Lagrangean,  are  considered.  In  this  process  two  new 
consistent  rate  complementary  energy  principles  for  the  rate  problem  of  finite 
strain  classical  elasto-plasticity  have  been  found.  Systematic  procedures 
to  exploit  these  complementary  rate  principles  in  the  context  of  assumed  stress 
finite  element  methods  are  also  discussed.  In  addition,  a critical  evaluations 
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of  the  relative  effectiveness  of  general  rate  principles,  in  both  TL  and 
UL  rate  forms,  and  in  terms  of  alternate  stress-rate  and  conjugate  strain-rate 
measures,  for  application  to  numerical  analysis  of  finite  strain  elasto- 
plasticity  problems,  is  made.  Also  included  in  the  present  paper  are  certain 
general  modified  variational  theorems  which  are  of  significance  in  the 
numerical  treatment  of  near  incompressible  behaviour  at  large  magnitudes  of 
plastic  strain. 

II.  Preliminaries: 

For  simplicity  we  refer  all  configurations  of  the  body  to  a fixed  rect- 
angular cartesian  coordinate  system.  We  adopt  the  notation:  (-)  under  symbol 
denotes  a vector;  (~)  under  symbol  denotes  a second-order  tensor;  a = A . b 

— M 

implies  that  a^  = A^  b^;  A . jB  denotes  product  of  two  tensors  such  that 

& • S>ij  = Aik  Bkj;  (£  : & = trace  (&T  •*£)  = AijBij;  and  0.  • I = Vi* 

A particle  in  the  undeformed  body  has  a position  vector  x = (x  e ) 

~ ono' 

where  e^  are  unit  cartesian  bases.  We  adopt  the  notation  v°  = d/kx^), 

in  the  undeformed  configuration,  Cq.  The  position  vector  of  the  particle 
in  the  current  deformed  state,  say  C^,  is  £ = (y^e^) . We  also  use  the  no- 
tation that  = (e^  3/dy^)  . The  gradient  of  is  the  tensor  F,  ie.,  F ■ 
o T 

(V  y ) ; F . =y.  = 3y./5x  . We  also  note  that  the  base  vectors  of  the 

~ i,o  r,o'  i a 

convected  coordinates  x^  in  the  current  deformed  configuration,  C^,  are  given  by, 

= .e^  a . The  non-singular  tensor  £ is  considered  to  have  the  polar- 

decomposition,  F = or  . (I  + h) , where  (I  + h)  is  a symmetric  positive  definite 

tensor  called  the  stretch  tensor,  ^ is  an  identity  tensor,  and  a is  an  ortho- 

T -1 

gonal  rotation  tensor  such  that  a - a . The  deformation  tensor  G is 
T 2 

defined  by  G 1 F . F = (I  + h)  . The  Green-Lagrange  strain  tensor  is 

defined  by  g = 1/2  (G  - I)  = 1/2  (e  + e^  + eT  . e)  where  e is  the  gradient 

o T 

of  the  displacement  vector  u (s  i - x) , ie.,  e = (V  u)  such  that  e:  = u. 

“ — ^ — lei'  i 
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For  our  present  purposes,  we  introduce  the  stress  measures:  (i) 

— ck'B 

the  true  (Cauchy)  stress  tensor  T (=  e^e ^ = T ; (ii)  a weighted 

stress  tensor,  generally  referred  to  as  the  Kirchhoff  stress  tensor,  a - 

cyB  _ cyB  r 

(=  a = w*lere  J i-s  determinent  of  the  matrix  L y^  ^J); 

(iii)  the  Piola-Lagrange  or  the  First  Piola-Kirchhof f stress  tensor,  ^t  = 

(t  .e  e.);  and  (iv)  the  second  Piola-Kirchhoff  stress  tensor,  s (= s ae  eQ) 
aj-rr-j  ~ cyP- cr"P 

As  discussed,  for  instance,  by  TRUESDELL  and  NOLL  (1965)  and  FRAEIJS  DE 
VEUBEKE  (1972),  the  relations  between  the  above  stress  measures  are  seen  to  be: 


-It,  --It, 

T T F . t _ F . 

^ J M M J M 


. F 


t = s . FT  = J(F_1  . t) 

M M M Vw  M 


(1) 

(2) 


and 


-1  -T  -T 

s = J(F  . T . F ) - t . F 

M ^ MM 


(3) 


-T  '-1  T 

where  F = (F  ) • We  also  note  that  eventhough  a and  ^s  are  distinctly 

cyB 

different  tensors,  they  have  the  interesting  property  that  a = s Q where 

OfP 

cyB 

the  contravariant  components  a and  the  components  s Q are  as  defined 

cyP 

earlier.  Finally  we  introduce  a stress  measure  which  is  labelled  by  FRAEIJS 
DE  VEUBEKE  (1972)  as  the  Jaumann  Stress,  r,  through  the  relation, 


r = 1/2  (t  . o'  + cyT  . tT) 

M M M M M 


= 1/2  s . (I  + h)  + (I  + h)  . s I 


W 

(5) 


It  is  seen  from  the  above  that  the  stress  tensors  £,  a,  £,  and  £ are  symmetric, 
while  t is  unsymmetric. 

M 

Also,  as  shown,  for  instance,  by  TRUESDELL  and  NOLL  (1965)  and  FRAEIJS 
DE  VEUBEKE  (1972),  the  linear  momentum  balance  (LMB)  equation,  the  angular 
momentum  balance  (AMB)  equation,  the  traction  boundary  condition  (TBC),  and 
the  displacement  boundary  condition  (DBC)  can  be  written  as  follows: 
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LMB: 


£N 

, N_ 

. T + p B = o 
r — 

(6) 

or 

7° 

• <£  • IT)  + P°B  - o 

(7) 

cr 

7° 

. t + p°B  • o 

N r — 

(8) 

N o 

where  B is  the  body  force  per  unit  mass,  and  p and  p are  the  mass  densities 
per  unit  volume  in  and  CQ,  respectively. 

AMB : 

T - tT  (9)  or  4 - / (10) 

or  F . t « tT  . FT  or  (h  + I)  . t . a * symmetric  (11a  ,b) 

f**  **  IN  N AS  A/ 

TBC:  "t  * n .'  t = n . (s  '.  FT)  ' * (12a ,b) 

' — 1 ^ »v  /v 

where  n is  a unit  normal  to  the  surface  Sa  , in  configuration  CQ,  where 

o 

tractions  are  prescribed  to  be  t per  unit  area. 

DBC ; 

u “ ti  on  S (12c) 

o 

where  S is  the  surface  in  C where  displacements  are  prescribed  to  be  u- 
u o 

o 

it  is  possible  to  have  a mixed-mixed  problem  wherein,  at  any  point  on  the 

surface  of  the  body  in  C certain  components  of  tractions  and  the  conjugate 

o 

components  of  displacements  may  be  simultaneously  prescribed. 

In  connection  with  the  rate  formulations  of  classical  elasto-plasticity,  the 
requirements  for  a suitable  stress  rate  are  now  well  recognized  as:  that  the 
stress  rate  vanishes  when  the  solid  continuum  undergoes  a rigid  motion  alone 
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and  whan  the  stress  tensor  is  referred  to  a coordinate  system  undergoing 
the  same  motion;  and  that  the  rate  of  invariants  of  the  stress  tensor  is 
stationary  when  the  stress  rate  vanishes.  The  questions  of  objective  stress 
rates  and  their  use  in  classical  rate  theories  of  plasticity  have  been  dis- 
cussed by  several  authors;  for  instance,  OLDROYD  (1950,  1958),  TRUESDELL  (1955), 
COTTER  and  RIVLIN  (1955),  PRAGER  (1961),  SEDOV  (1966),  MASUR  (1961), 

NAGHDI  and  WAINWRIGHT  (1961),  and  HILL  (1967).  With  this  background,  we  now 
discuss  the  following  rate  variational  principles  in  the  rate  theory  of 
classical  elasto-plasticity . 

III.  Finite  Strain  Elastic-Plastic  Analysis: 

3.1  Rate  Variational  Principles  in  Updated  Lagrangean  (UL)  Formulation 

In  the  UL  formulation,  we  refer  the  solution  variables  (displacements, 

strains,  and  stresses)  in  the  state  to  the  configuration  of  the  body  in 

the  immediately  preceding  state,  Cjj,  which  is  pressumed  to  be  known.  Let 
N 

yi  be  the  current  (in  state  C^)  Cartesian  spatial  coordinates  of  a particle, 

to  be  used  as  a reference  system  for  the  current  increment,  ie.,  from  to 
N 

Let  T be  the  true  stress  in  C„.  Thus  in  the  UL  rate  formulation,  one 
N+l  ~ N 

is  concerned  essentially  with  an  initial  stress  problem,  whereas  the  initial 
displacements  in  CN  as  referred  to  CN  itself  are  clearly  zero.  Let  ji, 
and  r represent  the  rates  of  the  Piola-Lagrange  stress,  2n<*  Piola-Kirchhof f 
stress,  and  Jaumann  stress,  respectively,  referred  to  the  current  configuration. 
Further,  we  note  that  s = s„  - t ; t = t_,  - T , etc.,  where  s„  is 

defined  as  the  2nt*  Piola-Kirchhoff  stress  in  state  as  referred  to  the 

configuration  (ie.,  measured  per  unit  area  in  CN) , etc.,.  Let  V represent 
the  gradient  operator  in  the  current  coordinates,  and  jj  be  the  rate  of  deforma- 
tion from  the  current  configuration.  We  define  the  rate  of  displacement 

gradient  e = (VNti)T  and  write  e = e + d),  where  did  = 1/2  (du./dy1?  + du./Sy1?)  I 

*■'**'•*'  **  l i j i j j l J 
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is  the  UL  Strain  rate  and  d>  = 1/2  (3u^/5y^  - Suj/dy^)J  is  the  spin-rate. 

3:1:0:  Rate  Potentials; 

★ 

We  also  use  d to  denote  the  corotational  rate  (or  what  is  also  usually 
called  the  "Jaumann  rate")  of  Kirchhoff  stress,  o.  Based  on  general  dis- 
cussions concerning  stress  rates,  contained,  for  instance,  in  the  References 


cited  earlier,  it  is  seen 

that; 

s = d - 6 

. tN  - 

tN 

. & 

(13) 

/■w 

/■w 

** 

and  £ = d*  - e 

r*j 

. tN  - 

tN 

. (!) 

(14) 

Considering  a classical  elasto-plastic  theory,  it  has  been  noted  by  HILL 
( 1967a, b)  that  a rate  potential,  V,  can  be  written  for  d*  such  that, 

d*  = dV/de  (15) 


As  also  noted  by  HILL  (1967a),  the  form  of  the  rate  potential  V can  be 
written  as : 


* - 1/2  ' 5 <Xk/k/  <16> 

which  yields  a bi-linear  constitutive  law  through  Eq.  (15).  Following  HILL 
(1967a)  we  note  that,  in  Eq.  (16),  L^^  Is  a tensor  of  instantaneous  elastic 
moduli! , assumed  to  be  positive  definite  and  symmetric  under  ij*— *kX  inter- 


change, a “ 1,  or  0 according  as  Is  positive  or  negative;  is  a 

tensor  normal  to  the  hyperplane  interface  between  elastic  and  plastic 
domain  in  the  strain-rate  space;  while  g is  a scalar  related  to  the  measure 
of  rate  of  hardening  due  to  plastic  deformation.  Prandtl-Reuss  type  rate 
equations  of  the  form  of  Eq.  (15)  for  classical  isotropically  hardening 
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materials  have  been  used  by  several  authors,  for  instance,  HUTCHINSON  (1973), 
McMEEKING  and  RICE  (1C  7 5)  , and  NEMAT -NASSER  and  TAYA  (1976). 

From  Eqs . (13)  and  (15)  it  is  seen  that  if  a rate  potential  V exists 

* • • 
for  d , then  a potentialW  also  exists  for  s,,  such  that. 


s = 3W/3e 


and  further  it  is  seen  that 


W = V - T : (£  . £) 


Likewise,  from  Eqs.  (14,  15  and  13),  it  can  also  be  seen  that  a rate  potential 
U for  t also  exists  such  that, 

£ = 3tf/eT  (19) 

where  U = V - t^:  (e  . £)  + 1/2  T^:  (e^  . e)  (20) 

= W + 1/2  tN:  (eT  . e)  (21) 

Further,  by  applying  the  polar-decomposition  theorem,  we  see  that. 


N+i 

£n 


_ /V7n  n+i.t  _ n+i  /T  . 

= (.11  > = £n  • + 


Where  F^!+^  and  <x!+^  are  the  deformation  gradient  and  rotation  tensors,  respec- 
~N  ~N 

tively,  in  C.,,.  as  referred  to  CM;  and  h is  the  UL  rate  of  stretch.  Writing 
N+l  N ** 

N+l 

= ^1  + g (where  a is  the  UL  rate  of  rotation),  it  is  seen  from  Eq.  (22) 

that, 

e = a + fi  (23) 

Thus,  as  may  be  expected,  it  is  seen  that  in  the  UL  rate  formulation,  fi  = & 
and  a = cl)  where  6 and  il)  are  as  defined  earlier.  From  the  definition  of  the 

^ i>w  rw 

Jaumann  stress  r,  as  given  in  Eqs.  (4,5),  the  UL  rate  i is  seen  to  be  given  by; 
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i - 1/2  + iT  + l • z + z • zNi  <2Aa> 

H 1/2  [£  + £T  + XN  . cfi  + d,T  . tN]  (24b) 

or,  also,  r = s + 1/2  (tN  . i + i,  . tN)  (25) 

comparing  Eqs . (13)  and  (25),  it  is  seen  that, 

r - i*  - 1/2(4  . t,"  + tN  . i)  (26) 


* ^ 

Thus  if  V is  a rate  potential  for  a , it  then  follows  from  Eq.  (26)  that 
there  exists  a rate  potential  Q for  r,  such  that 


r = 3Q/3H  = dQ  fbi 

(27) 

where , 

Q = V - 1/2  tN:  (£  . €) 

(28) 

Results  given  in  Eqs.  (27)  and  (28)  are  useful,  as  shown  later  on,  in  formulat- 
ing consistent  complementary  energy  rate  principles  for  the  rate  theory  of 
finite  strain  plasticity. 

With  the  above  considerations  for  rate  potentials  for  alternate  stress- 
rates,  we  proceed  to  formulate  the  boundary  value  problem  in  terms  of  the 
piecewise  linear  incremental  field  equations  and  boundary  conditions  in  the 
UL  rate  form,  as  follows: 

3.1.1:  Field  Equations  and  Boundary  conditions  in  UL  Rate  Form; 

• • • 

In  terms  of  s , e and  u 

^ /v  — 

Considering,  for  instance,  the  linear  momentum  balance  equation  for 

N+l 

s„  in  the  UL  Coordinates,  we  see  from  Eq.  (7)  that  this  can  be  written  as: 


_,n  r n+i 

7 . r 


L~N 


<£n+1)T]  + e*  £N+1  = 0 


(29) 


where,  as  defined  earlier,  is  the  2n<*  Piola-Kirchhof f stress  in 

N+l  N " 

as  measured  per  unit  area  in  C„  and  s.,  = T +s  etc.  The  initial  stress 

N 

field  £ , which  is  assumed  to  be  equilibrated  (in  an  actual  numerical  im- 
plementation this  may  not  be  true,  and  hence  it  is  often  necessary  to  employ 
corrective  iteration  procedures  to  check  the  true  equilibrium  in  C„),  is 
then  required  to  satisfy  the  equation. 


N N N N 
V . T + p B = o 

— r+t  r — 


Comparing  Eqs . (29) iand  (30)  we  obtain  the  rate  form  of  the  linear  momentum 
balance  equation  for  £.  Using  arguments  similar  to  the  above,  the  following 
field  equations  and  boundary  conditions  are  derived  in  UL  rate  form. 


(IMB)- 


(AMB)- 


v*  . [s  + tN  . (VNu3  + PNB  = 0 


* • — 
_ T 
s = s 


(Compatibility)-*  J = & = 1/2 (je  + eT)  = 1/2^  (V^u)  + (VNji)TJ 
(TBC)~*  n . [s  + t . (V  u)]  s t = t at  S 

— “ **  t s.  ■—  — — — Q 


(DBC)- 


u = u at  S 


N * 

where  p is  the  mass  density  in  Cn;  B are  rate  of  body  forces/unit  mass; 

, and  are  appropriate  segments  of  the  boundary  cf  the  solid  in  CN; 
and  n is  a unit  normal  to  the  boundary  of  the  solid  in  C^. 

In  terms  of  £,  e,  and  u 

In  manner  analogous  to  above,  these  field  equations  can  be  shown  to  be: 


(LMB)-* 


(AMB)- 


N , . N* 

V . t + p B = o 

(VNu')T  . TN  + £ = tT  + tN  . (A) 


or,  equivalently,  (AMB)  can  also  be  written  as: 


(AMB) 


. tN 

+ h . 

tN 

+ t 

= tT 

N 

+ T 

. h 

+ tN 

** 

N 

N 

• 

• 

*T 

N 

• 

•T 

. T 

+ T 

• SL 

+ t 

= t 

— 

+ T 

• h 
** 

+ z 

(Compatibility)-*  = (V*u_)'r 

or,  equivalently,  ( u_) T = a + h 


(TBC)-* 

ic  JL 

n . t = t at  S 

— M f"T 

r-w  U 

N 

(DBC)- 

same  as  in  Eq.  (35). 

(38) 
(38a) 

(39) 

(40) 

(41) 


3.1.2  General  Variational  Principles  in  UL  Rate  Form: 

3. 1.2.1:  In  terms  of  s,  £ and  u: 

Using  a virtual  work  principle  as  applicable  to  an  initial  stress  pro- 
blem, and  following  the  procedure  outlined  by  WASHIZU  (1975),  we  obtain  a 
general  UL  rate  principle  for  elastic-plastic  problems,  analogous  to  the 

well  known  HU  (1955)  - WASHIZU  (1955)  principle  of  linear  elasticity.  This 

* 

general  rate  principle  governing  Eqs.  (31  to  35  and  17)  can  be  stated  as 
the  condition  of  stationarity  of  the  functional: 

= J {w<£>  - PNB  • u + 1/2  tN:  [(VNu)  . (A)T] 

VN 

- s:  [e  - 1/2  < (v\)  + (^U)T  >]}dv  - J t . u ds 

“ J .t  . (u  “ u)ds . (42) 


where  W is  a rate  potential  for  £,  as  defined  through  Eqs.  (18  and  16).  The 
above  rate  variational  principle  governs  the  rate  variables  from  CN  to 

*The  general  variational  principles  as  stated  in  Eq.  (42)  above,  as  well  as 
those  in  Eqs.  (53,  69,  82,  85,  114,  119,  and  127)  to  follow,  can  be  modified 
appropriately  through  the  method  of  Lagrange  Multipliers,  as  discussed  in  ATLURI 
and  MURAKAWA  (1977),  to  account  for  discontinuities  at  interelement  boundaries 
when  these  principles  are  applied  to  a finite  element  assembly. 
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The  satisfaction  of  the  fully  nonlinear  field  equations  in  C^,  in  a numerical 
solution  method  such  as  the  finite  element  method,  must  be  checked  at  each 
step;  and  these  checks  can  be  performed  based  on  a variational  principle 
governing  the  nonlinear  field  equations  at  C^.  The  details  of  such  corrective 
iterative  procedures,  at  the  end  of  each  increment,  as  the  "equilibrium  check" 
"compatibility  mismatch  check",  etc.  can  be  found,  for  instance,  in  the  thesis 
by  MURAKAWA  (1978). 

We  now  consider  certain  special  cases  of  the  general  rate  principle 
given  through  Eq.  (42).  If  Eqs . (17,  33,  and  35)  are  met  a priori,  one  can 
reduce  Eq.  (42)  to  a functional  associated  with  the  UL  rate  principle  of 
potential  energy,  as: 

%2(i)  " J ’ {w(u)  - pNB  . u + 1/2  tN:  l(?u)  . (VNu)T]}dv 
P v ^ ■> 

N 

« 

- J t . uds  (43) 


The  principle  6np  = 0 can  be  seen  to  lead  to  Eqs.  (31,32,  and  34),  and  is 
equivalent  to  a principle  stated  originally  by  HILL  (1959). 

By  inverting  the  relation  in  Eq.  (17)  to  express  £ interms  of  £,  one 
can  achieve  a contact  transformation, 


(44) 


Using  Eq.  (44)  one  can  eliminate  £ from  Eq.  (42)  to  derive  a HELLINGER  (1914) 
REISSNER  (1950)  type  UL  rate  principle,  with  the  associated  functional: 


a J {-  £*<£>  - PN»  • u + 1/2  XN:  [(V*^)  . (VNi,)T] 
VN 

+ 1/2  £:  [(VNu)  + (VNu)T]|dv  - J _t  . _u  ds  - _t 


(ll  - jj)ds 


(45) 
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Based  on  the  arguments  presented  earlier  for  a linear  elastic  case, 
for  instance,  by  ATLURI  (1975),  it  can  be  seen  that  in  a finite  element 
application  of  the  principle  stated  through  Eq.  (45),  one  needs  to  assume 
over  each  finite  element,  an  arbitrary,  symmetric,  and  differentiable  stress- 
rate  field  £,  and  a differentiable  u that  is  also  inherently  compatible  at 
the  interelement  boundaries . 

We  now  examine  the  possibility  of  deriving  a complementary  energy  rate 
principle  from  Eq.  (45).  To  this  end,  we  first  note  that  in  the  present 
UL  rate  formulation  the  DIB  conditions  Eq.  (31),  are  linear  in  s and 
further,  unlike  in  Eq.  (7),  do  not  involve  coupling  of  £ with  displacement 
gradients.  Thus  it  becomes  possible  to  satisfy  both  the  1MB  and  AMB  conditions, 
Eqs . (31  and  32)  respectively,  a priori,  by  choosing  a symmetric  s such  that, 

s = curl  curl  A + sF  (46) 


where  A is  the  symmetric  Maxwell-Morera- Beltrami  second-order  stress  function 
tensor  for  a general  three-dimensional  case.  In  Eq.  (46)  curl  A is  defined 
such  that  (curl  A)^  = e^^,  A4L-  (curl  curl  A),..  = e4 e,_  A 


• 0 

imn  jpq"mp,nq*  ijk 


ipk  jk,p 

the  alternating  tensor,  and  s p is  any  symmetric  particular  solution  such 
that. 


_N  .p  _ N*  N r N /rjr.n 
1 • £ = -p  B - V . [t  . (\Tu)J 


(47) 


one  simple  way  of  satisfying  Eq.  (47)  is  to  assume  particular  solutions  for 
the  direct  stresses  s^  ( no  sum  on  i;  i = 1,2,3)  only,  in  the  following  way. 


• P _ p r N . N . v s 

sii  JN  ^"p  Bi  ' (Tkj  °i; j ’k^dyi 


(48) 


(no  sum  on  i,  1 = 1,2,3) 


and 


•Sj  - 0 


(i  + J) 


(49) 
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In  Eq.  (48)  ( ; i)  indicates  8(  )/dy^.  However,  if  the  above  assumptions  are 

l 

used  in  an  assumed-stress  type  numerical  scheme  (discussed  further  below) 
the  question  of  completeness  of  the  chosen  stresses,  or;  in  other  words, 
the  effect,  on  that  numerical  solution,  of  the  lack  of  account  of  the  influence 
of  displacement  rates  on  the  chosen  shear-stress-rate  field,  as  in  Eq.  (49), 
remains  to  be  answered.  Such  effects  can  be  only  understood,  in  general, 
from  a detailed  mathematical  study  of  convergence  of  the  method,  which  study 
is  not  pursued  in  the  present  paper  and  remains  an  open  question.  Assuming 
that  the  satisfaction  of  the  LMB  and  AMB  conditions  in  the  manner  of  Eqs.  (48  and 
49)  is  "satisfactory",  and  further  if  the  TBC  condition  is  also  met  a priori, 
one  can  reduce  Eq.  (45)  to  a functional  associated  with  the  complementary 
energy  principle. 


*1  (H>i>  = J {*  i*  (i)  - 1/2  £N : [(VNu)  . (VNu)T]}dV 


+ J _t  . ii  ds 

S„ 


In  a finite  element  application,  VN  can  be  subdivided  into  M subdomains, 

VmN  = *»••*  each  with  a boundary  In  general  it  is  seen  that 

dV  „ = S + S + p „ where  S„  , S are,  respectively,  the  portions 
mN  mN  mN  mN 

of  dV  „ where  tractions  and  displacements  are  prescribed,  and  o is  that 
mN  ’ rmN 

portion  of  dv  which  is  common  to  an  adjacent  element  (interelement  boundary). 
mN  n. 

It  is  to  be  noted  that  in  the  finite  element  application  of  Eq.  (50),  the 
candidate  stresses  s^  should,  not  only,  satisfy  the  IXB  and  AMB  conditions, 

Eqs.  (31)  and  (32),  but  also  satisfy  the  interelement  traction  reciprocity 

condition^,  (_t)+  + (_t)"  = 0 (where  t^  is  defined  through  Eq.  (34))  at 
a priori.  One  can  introduce  this  interelement  condition  directly  as  a condition 
of  constraint  into  Eq.  (50),  in  order  to  preserve  a wide  choice  of  candidate 


The  superscripts  (+)  and  (-;  denote,  respectively,  the  two  sides  o.  p in 

the  limit  that  p „ is  approached. 
mN 
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stress-rates  s , as : 


(H-i-l)  'I  J {-S*(4>  - 1/2  tN:  [(vNu)  . (/u)Tl}<lV 


m V 


mN 


+ L I t . u ds  + ^ J I . 2-  dS 

m S m p r 

M * I 


(51) 


mN 


mN 


where  are  Lagrange-Multipliers  to  enforce  the  constraint  of  interelement 
traction  reciprocity,  and  these  can  be  seen  to  be  the  displacement  rates  at 
the  interelement  boundaries.  The  basic  idea  of  choosing^s,  as  in  Eqs.  (47  and 
48),  to  satisfy  the  121B  condition,  and  the  modified  complementary  energy 
functional  of  Eq.  (51),  were  used,  in  a somewhat  less  general  fashion  than  is 
given  here,  by  ATLURI  (1973)  in  formulating  a finite  element  method  and 
applied  to  solve  the  problem  of  buckling  of  a shallow  arch.  The  short- 
comings of  this  approach  were  later  discussed  by  ATLURI  and  MURAKAWA  (1977). 

Moreover,  in  the  modified  complementary  energy  rate  principle  of  Eq.  (51), 

• • 

in  addition  to  the  fields  u and  £ within  each  element,  the  displacement-rate 
field  at  the  interelement  boundaries  also  enters  as  an  independent  variable. 
However  in  the  finite  element  application  of  the  Hellinger-Reissner  type 
principle  as  in  Eq.  (45),  only  the  two  variables  s (arbitrary,  but  differenti- 
able and  symmetric  second  order  tensor)  and  u (differentiable  and  interelement- 
compatible)  need  to  be  assumed  within  each  element;  and  the  interelement 
traction  reciprocity  condition  then  follows  a posteriori  from  the  variational 
principle.  Considering  this,  and  the  fact  that,  in  addition,  one  must  study 
the  convergence  properties  of  the  finite  element  scheme  based  on  Eq.  (51) 
to  assess  the  effects  of  choosing  the  particular  solution  in  a specific 
way  as  in  Eqs.  (48,  49),  it  appears  most  consistent  and  rational  to  directly 
apply  the  Hellinger-Reissner  type  principle  as  in  Eq.  (45).  The  finite 
element  method  thus  generated  based  on  assumed  s as  well  as  u,  can  be  called 
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a "mixed-method",  which  leads  to  simultaneous  algebraic  equations  for 

finite  element  nodal  displacements  as  well  as  nodal  stresses  (or  alternatively, 

the  nodal  values  of  the  stress  functions  A as  defined  in  Eq.  (46)). 

Finally,  it  is  noted  that  the  IiiB  condition  may  be  satisfied  more 
easily  by  choosing  £,  unlike  in  Eqs.  (46,47),  such  that, 

s - curl  curl  A - . (V^u)  + & (52) 

where  s*5  is  such  that  v . s^  ■ -p  B.  However  the  chosen  s as  in  Eq.  (52) 
then  ceases  to  be  symmetric,  and  thus  the  AMB  condition  must  be  introduced 
as  a constraint  condition,  into  the  associated  complementary  energy  functional 
of  the  type  given  in  Eq.  (50),  through  additional  Lagrange  Multipliers. 

Thus  it  appears  that  a rate  complementary  energy  principle  in  UL  form 
based  on  j may  not  be  consistent  and  practically  useful  in  the  analysis  of 
finite  strain  plasticity  problems. 

3. 1.2. 2:  In  terms  of  £,  e,  and  u 

Analogous  to  the  way  discussed  earlier,  it  can  be  shown  that  Eqs.  (36, 
37,39,41,35,  and  19)  follow  as  the  Euler  equations  and  natural  boundary 
conditions  corresponding  to  the  stationarity  of  the  functional. 


= J {u(e)  - PNB  . * + ,tT:  lA)T  - e]}dV 

VN 

• • • • • 

- J _t  .ju  ds  - J ji  . (u  - u)  ds 


where  U is  the  rate  potential  for  as  defined  through  Eqs.  (20,21). 

We  now  consider  certain  special  cases  of  the  above  principle.  If  Eqs. 

• • 

(19,39,  and  35)  are  met  a priori,  one  can  eliminate  £ and  JC  as  variables  from 
Eq.  (53)  and  derive  a rate  functional  governing  the  rate  potential  energy 
principle,  as: 
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(54) 


• C 9 * • • 

tt  (u)  « J |u(u)  - pNB  . ujdV  - J*  t . uds 


N 


'N 


This  rate  variational  principle  was  first  stated  by  HILL  (1959)  and  has 
been  widely  used  in  finite  element  applications  to  elastic-plastic  problems 
[See  for  instance,  NEEDLEMAN  (1972) , McMEEKING  and  RICE  (1975),  and  NEMAT- 
NASSER  and  TAYA  (1976)]. 

It  is  interesting  to  note  that  both  the  IMB  and  AMB  conditions, 

Eqs.  (36)  and  (37)  respectively,  as  well  as  the  TBC,  Eq.  (41),  must  follow 

★2  * *2 

from  the  principle,  Stt^  (6u)  = 0,  with  tt^  given  as  in  Eq.  (54).  It  is 
shown  below  that  the  AMB  condition  is  inherently  embedded  in  the  special 
structure  for  U.  Thus,  using  the  definition  of  t (=3u/3er)  and  Eq.  (21), 
it  is  seen  that. 


oe  oe 


&W  N *1 

= + t . e 

de  3eT  ~ 


(55) 


(56a) 


. N T 
s + T . e 


(56b) 


Wherein,  the  definition  of  s from  Eq.  (17)  has  been  used.  Substituting  for 

^t  from  Eq.  (56b)  into  the  AMB  condition,  Eq.  (37),  it  is  seen  that  the  AMB 

• 

condition  is  inherently  met.  This  is  due  to  the  special  structure  for  U 

as  given  through  Eq.  (21).  Conversely,  it  is  seen  that  if,  instead  of 

• • 

Eqs.  (20,  and  21),  an  arbitrary  U is  postulated  as  a function  of  £,  then  the 
principle  based  on  the  functional  in  Eq.  (54)  ceases  to  be  valid,  since  the 
AMB  condition  ceases  either  to  be  built  into  the  structure  of  U or  to  follow 
unambiguously  as  an  Euler  equation  from  the  vanishing  of  the  first  variation 
of  the  said  functional.  This  fact  appears  to  be  never  explicitly  stated  in 
the  literature. 
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Now,  by  inverting  the  bi-linear  relation  in  Eq.  (19)  one  may  achieve  a con- 
tact transformation. 


U 


.*  . 
-E  (t) 


(57) 


such  that,  — = e*  (58) 

hi  ~ 

Using  Eq.  (57)  to  eliminate  e from  Eq.  (53)  one  may  formally  obtain  a 
functional : 

nHR  <£»“>  " J*  {-  E*(t)  - pj)  . u + tT:  [(VNu)T]}dV 
VN 

• • • • • 

- J*  _t  . u ds  - J £ . (u  - u)  ds.  (59) 


In  fact  the  above  functional,  as  the  basis  for  a Hellinger-Reissner  type 

variational  principle  was  used  by  NEALE  (1972).  However,  the  validity  of  such 

a principle  needs  a closer  examination.  If  &tt  = o,  with  rr  as  given  in 

Eq.  (59),  is  a valid  Hellinger-Reissner  type  rate  principle,  we  note  that  the 

corresponding  Euler  equations  and  natural  boundary  conditions  must  be:  (i) 

the  LMB  condition,  Eq.  (36);  (ii)  the  AMB  condition,  Eq.  (37);  (iii)  the 

compatibility  condition,  Eq.  (39),  (iv)  TBC , Eq.  (41);  and  (v)  the  DBC , Eq.  (35). 

It  is  seen  upon  examining  Eqs.  (17)  to  (25)  of  NEALE'S  (1972)  development, 

the  AMB  condition,  Eq.  (37)  of  the  present  paper,  is  in  fact  not  an  Euler 

equation  of  the  principle  61  = 0 [Eq.  (22)  of  NEALE'S  (1972)  work,  which  is 

identical  to  6tt  =^0  of  the  present  article].  Thus,  if  at  all  the  AMB 

condition  is  satisfied,  for  the  validity  of  the  principle,  this  condition  must 

be  embedded  in  the  special  structure,  if  any,  for  the  complementary  energy 
**  * 

density  function,  E (t)  of  Eq.  (57)  above.  To  examine  this  possibility, 


20 


consider  the  form  of  U as  given  from  Eqs.  (20  and  16): 

2U  = «ijLijkiekjt  -2  | -2  Tijeik*kj  + TijUk,iUk,j  (60) 

The  stress  rate,  t,  as  derived  from  the  above,  is 


„ 3U  _ „ a > ^ kX  eXk^  , V_N 

2t. , = 2 — ; (L.  „ -2  — X X.  .) . - (e.  + e . )t  . 

ij  g ijkX  g ij  kX  2 v im  mi7  mj 


ji 


N , * 

‘ Tix(exj  " ejx* 


(61) 


where  I^jkjj*  O',  g and  \t#  are  as  defined  before.  The  constitutive  law,  Eq. 


kX 

• • 
(61),  is  of  bi-linear  type  for  t^  in  terms  of  e^* 


The  inversion  of  this 


relation  in  closed  form,  to  express  e^  in  terms  of  t^j  appears  to  be  impos- 
sible, in  general. 

Rearranging  the  right  hand  side  of  Eq.  (61),  it  can  be  rewritten  as: 


* 

j i ^ j ikn  ekn 


"k 

ot  matrix  form  {^}  9x1  t L^9x9  ^9x1 


(62) 

(62a) 


Where  dePends  on  the  tensor  of  instantaneous  elastic  modulii 

(which  has  the  symmetry  properties  = ^jikX  = ^jiXk  = ^kiij^  and  ot^er 

N ★ 

quantities,  a,  g,  and  t^.  However  Lj^n  has  the  only  symmetry  property. 


k * 

Ljikn  ~ Lknji 


(63) 


Eventhough  an  analytical  form  for  the  inverse  of  the  matrix  [ L]  of  Eq.  (62a) 
appears  impossible  to  be  obtained,  one  may  numerically  invert+  Eq.  (62,  62a) 
to  write. 


-fflowever,  in  the  "first  step"  of  the  solution  if  the  solid  remains  elastic, 
^nd  the  initial  stresses  £ are  zero,  in  the  initial  configuration  C , then 
L.  , 2 L.,..  and  *lence  9x9  matrix  of  Eq.  (62a)  cannot  be  inver?ed,  be- 
cause of  tne  symmetry  properties  of  stated  above. 
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* * -1  * r*  1-1  r'i 

eij  = Lijk4  tkX  or  [S)  = [ L]  {^} 


(64) 


* -1  * -1 

where,  in  general,  L^kjj  = LkXij*  UsinS  Ecl*  (64)  a contact  transformation 

**  * 

can  in  fact  be  made,  to  find  E (t)  such  that 


3E  /dt^  e Lijki  tkX 


(65) 


If  the  AMB  condition,  Eq.  (37),  is  inherently  embedded  in  the  structure  of 

•*  * 

E (^t) , then  the  condition. 


N 

e.„  t + t._  = symmetric 
ij  jk  ik 


(66) 


must  be  identically  satisfied  when  e..  is  expressed  in  terms  of  t through 

i j mn 

Eq.  (65).  Doing  so  we  see  that  the  AMB  condition  is  expressed  by  the  necessary 

**  * 

condition  on  the  structure  of  E (^t)  that 


or 


•*  ■ • 

(BE  ^dt^)  T^k  + t^k  must  be  symmetric 

* -1  n 

L.  . t T..  + t.,  must  be  symmetric 

ljmn  mn  jk  ik  1 


(67) 


1 II 

It  can  be  seen  that  neither  of  the  two  terms,  I and  II  above,  is  by  itself 

symmetric  under  i *—•  k interchange.  The  other  possible  ways  in  which  I + II 

above  can  be  symmetric  under  i *— • k interchange  are:  (i)  firstly,  one  term 

is  a transpose  of  the  other;  however,  it  is  easy  to  see  that  this  is  not  the 

case;  (b)  secondly,  the  first  term  can  be  expressed  as  the  sum  of  a symmetric 

* - 1 

tensor  and  the  transpose  of  the  second  term.  Since  L^k^  cannot  be  found 

analytically,  and  with  the  only  knowledge  that  = ^kXij’  aPPears 

impossible  to  verify  this  assertion. 

Thus,  since  the  AMB  condition  is  neither  clearly  an  Euler  equation 

*2 

corresponding  to  the  stationarity  of  tt  of  Eq.  (59),  nor  can  be  verified 

HR 

_ ** 

to  be  embedded  in  the  structure  of  E (t),  the  Hellinger-Reissner  type  principle 
based  on  Eq.  (59)  appears  to  be  of  little  practical  value.  For  similar 
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reasons,  the  complementary  energy  rate  principle  as  stated  by  HILL  (1959) 

(which  can  be  derived  formally  from  Eq.  (59)  by  requiring  t to  satisfy 

the  LMB  condition,  Eq.  (36),  and  the  TBC,  Eq.  (41),  a priori)  also  appears 

to  be  of  little  practical  value. 

• • • • • • • 

3. 1.2.3:  In  terms  of  r(t,or);  a^sco):  h,  u 

To  seek  alternative  ways  to  avoid  the  above  discussed  difficulties 

in  formluating  a consistent  complementary  energy  rate  principle  and 

Hellinger-Reissner  type  rate  principle,  we  transform  the  general  variational 

• • • • 

principle  associated  with  Eq.  (53)  into  one  involving  r,  a,  h and  u as 

A*  /V  /V 

variables. 

First,  by  comparing  Eqs.  (20)  and  (28)  we  note  that, 

U = Q + 1/2  tN:  (e.e)  - tN:  (e.e)  + 1/2  tN:  (eT.e) 


Q + 1/2  T»:  feTi)  + tN:  faT.e) 


(68) 


Using  (68)  to  express  U in  terms  of  Q (which  is  a function  of  h s e)  and 

A/' 

• • • 

writing  e = h + or,  we  rewrite  Eq.  (53)  as: 

***  *+  A# 

= 1 {QQ3>  + CarT.ar)  + rN:  (aT.h) 

VN 

- pNB.u  + tT:  [(VNu)T  - h - ff]}dV 

— — rw  A#  J 


^ • • • • 

- J*  t.u  ds  - J t^/u-u)ds . 


(69) 


N N 

In  the  above  or  is  a skew- symmetric  tensor.  The  condition  of  vanishing  of 
the  first  variation  of  the  above  functional  can  now  be  written  as: 
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6tt*2,(6u;  6h;  6a;  6t)  = f {[^  - l/2(t  + tN.  a + tT  + aT  . tN)]  : 6h 
HW  — ~ ~ ~ J ~ ~ ~ ~ ~ ' ~ 

V an 

N ~ 

+ C(VNu)T  - a - h]  : 6 tT  - [tN.  a f h . tN  + t]  : 6aT  - [vN.t  + pNBj  . 6uW 


i ^ t * « » 

- J ( t_  - n . ^t)  . 6u  ds  - J 6 . (u  - u)  ds  = 0 (70) 

S ~ S 

aN  UN 

Noting  that  a,  and  hence  6d,  are  skew-symmetric  tensors,  it  can  be 
clearly  seen  that  the  Euler  Equations  and  natural  boundary  conditions 
corresponding  to  Eq.  (70)  are:  (i)  the  constitutive  law,  Eq.  (27);  (ii) 
the  LMB  condition,  Eq.  (36);  (iii)  the  AMB  condition,  Eq.  (38);  (iv)  the 
compatibility  condition,  Eq.  (40);  (v)  the  TBC,  Eq.  (41);  and  (vi)  the 
DBC,  Eq.  (35). 

One  can  now  invert  (even  if  numerically)  the  relation  of  Eq.  (27)  and 
achieve  a contact  transformation, 


• • • * ic  * 

Q - r : h = -R  (r) 


• • •«!  N # *T  N •*  • 

or  Q - l/2[t  + t + t . or  + or  . t ] : h = - R (r) 

'V  /v  /v  /v  /v  «%/ 


such  that 


dR  for  = h 


(71a) 

(71b) 

(72) 


Substituting  Eq.  (71)  in  Eq.  (69)  we  can  derive  a functional,  involving 
• • • • 

only  u,  and ^r  (and  hence ^t  and  a)  as  variables,  corresponding  to  a Hellinger- 
Reissner  type  principle  which  has  as  its  Euler  Equations  and  natural  b.c, 

Eqs.  (36,  38,  40,  41,  and  35).  If,  in  addition,  one  assumes  that  the  LMB 
condition  and  TBC  for  t,  Eqs.  (36,  and  41)  are  satisfied  a priori,  one  can 
eliminate  h and  u as  variables  from  the  functional  in  Eq.  (69)  and  thus 
obtain  a complementary  energy  functional: 
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n£2 <&•’£>  = I {■  R*<i>  + 1/2  iN;  - i^jdv  + J (£*%$)  • £ 


* **T  N**T  m 

In  the  above,  the  definition  of  r [=  1/2  (t  + t + t . a + a1  . T )]  is 

implied;  and  the  spin-rate  field  s is  required  to  be  skew-symmetric.  The 

*2 

variational  equation,  6tt  = 0,  for  constrained  6t  (which  obey  the  constraint 

\j  ft* 

N * * * • 

V . 6t  = 0 in  V and  n . 6t  = 0 at  S ) and  for  constrained  6o  (which 
— ~ N - ~ o'  ~ 

N 

is  required  to  be  skew  symmetric)  is  seen  to  lead  to 


6tt^ (6^,60)  = 0 = J {[(VNu)T  - a - j:  6^tT  - [t  + h . tN  + *N  • «]  : 

d r 


SR*  1. 


“T'l  « * • i • 

a j-  dV  + J (n  . 6^t)  (u  - u) 

J r*  ' ~ 


Noting  that  by  definition,  dR  /d^r  = h,  it  is  clearly  seen  that  Eq.  (74)  leads, 

as  its  Euler  equations  and  natural  b.c,  (i)  the  compatibility  condition, 

Eq.  (40);  (ii)  the  AMB  condition,  Eq.  (38);  and  (iii)  the  DBC,  Eq.  (35). 

Thus,  Eq.  (73)  forms  the  basis  of  the  most  consistent  and  practically 

useful  rate  complementary  energy  theorem  for  the  UL  rate  formulation  of 

finite  strain  plasticity  analysis  methods  because:  (a)  the  admissible  t 

is  required  to  satisfy,  ^ priori,  only  the  uncoupled,  linear  LMB  equation, 

Eq.  (36),  and  TBC,  Eq.  (41),  which  can  be  met  easily  in  applications,  by 
' N 'p 

setting  t = V x Y + t where  Y are  first  order  stress  functions  (once- 

• N # 

differentiable)  and  t is  any  particular  solution  such  that  V . = 

N” 

-p  B;  (b)  the  AMB  conditions,  the  compatibility  condition,  and  the  DBC 
follow  unambiguously  as  Euler  equations. 

In  a finite  element  application  of  the  complementary  rate  principle 
as  stated  through  Eq.  (73),  the  assumed  stress-rate  field ^t  must  not  only 
satisfy  the  LMB  condition  (Eq  • 36)  within  each  element,  but  must  also 


satisfy  the  traction  reciprocity  relation  at  the  interelement  boundary 
* ’ + * * . - r 

viz.,  (n  .^t)  + (n  .^t)  at  [where  + and  respectively,  indicate  the 

two  sides  of  in  the  limit  that  pmM  is  approached].  This  may,  in  general. 


mN 


pose  a severe  restriction  on  the  choice  of  ^ within  each  element,  especially 
when  the  element  is  of  a arbitrary  curved  geometry.  In  such  a case,  it  may 
be  preferrable  to  include  this  interelement  traction  reciprocity  condition 
as  a constraint  condition  directly  into  the  functional  of  Eq.  (73).  The 
Lagrange  multipliers  introduced  to  this  end  can  be  seen  to  be  the  interele- 
ment boundary  displacements.  The  thus  modified  complementary  energy  rate 
principle  for  an  assembly  of  a finite  number  of  elements  can  be  stated 
as  the  stationarity  condition  of  the  functional: 

*2  ' ' ~ „ V fi-  r ’*  * , N *T  * ’T  ’1 

nMc  fe’-S’V  =lJ  L"  r <£)+ 1/2  z ■■  fe  •«)  - £ : sddv 


m V 


mN 


+ f (n  ,t)  . u ds  + f (n  .t)  . ^ ds^ 
— ~ — « — ~ nsJ 


(75) 


mN 


mN 


In  the  above  functional,  a and  ^t  are  chosen  independently  within  each 

element  in  terms  of  undetermined  parameters,  whereas  are  chosen  in  terms 

~P 

of  displacements  at  nodes  of  a finite  element  and  hence  u^  are  common  to 
elements  sharing  a common  boundary.  Thus  the  undetermined  parameters  in 
the  field  functions  for  o:  and  t can  be  eliminated  at  the  element  level  and 


expressed  in  terms  of  the  generalized  nodal  displacement  coordinates.  The 
finite  element  method  based  on  Eq.  (75)  thus,  in  the  end,  results  in  a 
standard  stiffness  matrix  procedure  [See  MURAKAWA  and  ATLURI  (1978a, b)  for 
instance,  for  details  of  finite  element  application  of  the  complemento,-y  energy 
rate  principles  in  finite  elasticity].  Alternatively,  the  interelement 
traction  reciprocity  can  be  satisfied  a priori  by  an  appropriate  choice  of 
the  first-order  stress  functions]?  from  which  the  equilibrated  t are  derived. 
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The  finite  element  method  then,  in  general,  will  lead  to  a "flexibility  matrix" 
type  approach.  Such  "stiffness"  and  "flexibility"  type  finite  element  methods 
based  on  Eqs.  (75)  and  (73),  respectively,  for  analyzing  certain  metal- 
forming problems  are  the  subjects  of  the  author's  work  in  progress  and  will 
be  subjects  of  a forthcoming  paper. 

3.1.3;  Near  Incompressibility  in  the  Fully  Plastic  Range: 

As  discussed  in  the  introduction,  an  important  aspect  of  numerical 

schemes  for  finite  strain  elastic-plastic  analysis  is  the  problem  of  accurate 

treatment  of  nearly  incompressible  deformation  rates  at  auch  magnitudes  of 

strain.  In  a finite  element  application  of  the  potential  energy  rate 

formulation  of  the  type  given  by  Eqs . (43  or  54) , if  the  assumed  deformation 

rates  do  not  a priori  obey  the  incompressibility  constraint,  it  may  be 

necessary  to  retain  this  constraint  as  an  a posteriori  constraint  through 

a Lagrange  multiplier,  the  hydrostatic  pressure.  To  this  end,  consider  the 
• • 

rate  potential  W(e ) , Eq.  (18),  for  a classical  Prandtl-Reuss  type  rate 
constitutive  law: 


W(e)  = 1/2  a 


(76) 


The  corotational  rate  of  Kirchhoff  stress,  a , for  a classical  Prandtl-Reuss 
type  approximation  can  be  written  in  terms  of  £ , as  suggested  by  McMEEKING 
and  RICE  (1975),  as 


•* 

a 


,N  ,N 

" = 2ur5  6 - — T_LLMi*  + x ’ 6 

ij  ZM,L  ik°j£  (2h  + 6 p.)  (-N  2 J®kX  X kk°ij 


(77) 


■ ' 'N 

where,  a * 1 if  at  yield  and  > 0 and  a = 0 otherwise;  is  the  devia- 

-N  2 

toric  Cauchy  stress  in  CN;  (T  ) * (3/2)t^t^;  h is  the  slope  of  the  uniaxial 

stress/plastic  strain  curve;  and  \ and  p,  are  Lame'*  constants.  We  rewrite 
Eq..  (77)  as: 
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- 1 


(78) 


where  the  definition  of  E^j^  *-s  apparent  from  comparing  Eqs.  (77)  and  (78). 
•* 

We  can  write  a in  terms  of  its  deviatoric  and  hydrostatic  parts,  as  : 
'*■*•*'  • 

a . . = a + (1/3)  a.  . 6 , . ; and  likewise,  express  e . . = e ! , + (l/3)e.  .6,  . 
ij  ij  kk  ij  r ij  ij  kk  ij 

From  Eq.  (78)  one  can  write; 

*ic 

ct.  = 2(1  E. 


O*  = (3 X + 2(X  )e 


ij  - ^ijlcTk X ’ kk 

Using  Eqs.  (77-79a,b)  in  Eq.  (76)  we  can  write. 


W + 2 <ekk)2  -l:  fe  • 


kk 


(79a, b) 


(80) 


or,  equivalently, 


W = ^EijWeWeIj  + 


3X  + 2 


^(ekk)2  - T^:  (c  . e) 


Vv  /v' 


(81) 


To  obtain  numerically  accurate  solutions  in  situations  of  fully  developed 

plasticity  it  may  be  advantageous  to  retain  the  hydrostatic  pressure  as  an 

independent  variable,  and  thus  derive  a mixed  variational  principle  which 

represents  a modification  to  the  potential  energy  rate  principle  given 

through  Eq.  (43),  with  W expressed  as  in  Eq.  (80).  To  this  end,  we  intro- 
• • 

duce  both  e„  and  as  additional  independent  variables  into  the  function- 
al of  Eq.  (43),  through  the  introduction  of  Lagrange  multipliers  |3  and  a, 

• • • 

respectively.  From  this  general  variational  principle  (with  u,  e^,  e^, 

• 

g , and  a,  all  as  variables)  we  demand  all  the  necessary  field  equations 
in  the  case  of  near  incompressibility.  When  the  Lagrange  multipliers  are 
identified  with  the  relevant  stress  rates,  this  general  rate  variational 
principle  can  be  written  as: 
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2*/  * * **  •*  v 

"riw^5  eij;  ekkJ  CTij*  CTkk> 


= £ Eijkiekjeeij  + 2(ekk)  + *-Uk,j6kj  ' ekk^ 

N 


k,  j kj  kk  (3X  + 2\x  ) 


* • • JL,  ACT.  1 0 , . • • 

» r ^ ir_*  kk  i j-,  N 

^U(i» j)  ’ ir*-'  ij  ' 3 \ + 1 Jriju(k,i)u(k,j) 

N * # N*  # T 

+ 1/2  T^u^  j “ P Biui  ] dV  + Boundary  terms 


The  above  general  principle  can  be  seen  to  be  valid  in  both  the  cases  of 


near  and  precise  incompressibility.  In  the  above,  the  notation,  u 


(k,i) 


1/2<"k,i+  Vk>-  and  v k = has  been  used.  If  from  Eq.  (82)  one 

• • • 

eliminates  (i)  by  defining  a priori  = u^  , and  (ii)  eliminate 


e.  . through  a contact  transformation: 


X * 2 • *ukk 

2 '■ekk;  " ekk  (3X  + 2 jjl) 


(3X  + 2 [i)‘ 


one  obtains  a mixed  variational  principle,  involving  u^  and  ct  as  variables. 


governed  by  the  functional. 


*2  / '*  ? f * * x ^kk' 

(ui;  CTkk}  "I  l^Ei  jki.u(k,X)u(i,  j)  " 2 ~ ~ 2^2 


/*  .2 
<CTkk> 


X •*  * N ' 

+ (3X  + 2n)  Gkku£,j6£j  " TijU(k,i)U(k,j) 


+ 2 TijUk,iUk,j  ‘ PViJdV  ■ S tiuids 

C 


k :.v 


which  remains  valid  for  nearly  or  even  precisely  incompressible  behavior  at 
large  plastic  strains  for  all  assumed  displacement  rates  u^  that  do  not  obey 
the  constraint  of  incompressibility  a priori.  Eq.  (84)  and  the  associated 
variational  principle  are  analogous  to  the  ones  in  the  case  of  linear 
isotropic  elasticity  given  directly  by  HERRMANN  (1965);  however,  the  way 
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in  which  HERRMANN  arrived  at  his  principle  for  the  linear  elastic  infinitesimal 

deformation  case  is  not  evident  from  reading  his  work. 

Likewise,  using  the  definition  of  W as  in  Eq.  (81)  in  Eq.  (43),  and 
• • 

introducing  and  as  additional  independent  variables  into  the  function- 
al in  Eq.  (43)  through  appropriate  Lagrange  Multipliers,  one  can  derive 
another  alternate  general  variational  principle,  which  remains  valid  in 
the  limit  of  incompressibility,  with  the  associated  functional: 


*2  .*  *,  • **,  **  . 

"hw  (v  'ir  ekk;  °ij!  “W 


•* 

a. 


* { Kju'Vi* + ^H^'kk + [“k,k  - «kiP  -r 

N 


r*.  t N ’ * IN'* 

+ l_u,,  - e : o . ' - T ,u.  . u.,  ,.+•=•  t . .u.  .u, 

(i.j)  ij  ij  ij  (k,i)  (k, j)  2 lj  k,i  k,j 


• • ) 

PNBiUiJ  ^ + boundary  termS 


(85) 


where  u^  ^ = u^  ^ - u^  ^6 ^ j / 3 and  u^  ^ is  defined  earlier.  If  from 

Eq.  (85)  one  eliminates  as  a variable  through  a priori  satisfying  the 
• • • 

condition  ejj  * u^  and  is  eliminated  through  the  contact  trans* 


formation  given  below. 


• •*  • 2 
.3X  + 2 u.  ‘ 2 €kk°kk 
1 6 ; ckk  " 3 “ 6(3X  + 2n) 


kk 


(86) 


one  obtains  an  .alternate  mixed  variational  principle,  also  involving  u^  and 


0.  as  variables,  governed  by  the  functional: 

Ivlv 

*2  / •* 


•* 

a 


VC  £ # . [ - 1,  It  * 

TTmP2(ui;akk)  = J l**EijkJtU(i, j)U(k,jt)  + T Uk,k 


N 


/*  x2 

■ (akk> 


6(3X  + 2m. ) ‘ TijU(k,i)U(k,j)  + 2 TijUk,iUk,j  " PnV,iJdV 


’I 


vVid8 

T 


(87) 
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Equation  (87)  and  the  associated  principle  are  analogous  to  those  derived 
by  KEY  (1969)  for  linear  elastic  infinitesmal  deformation  problems;  except 
that,  KEY  (1969),  derives  such  a principle  through  FRAEIJS  DE  VEUBEKE'S 
(1951)  interpretation  of  what  is  generally  known  as  the  Hellinger-Reissner 
theorem  in  linear  elasticity. 

NAGTEGAAL,  PARKS,  and  RICE  (1974,  Appendix  II  therein),  to  improve  the 
accuracy  of  UL  rate  finite  element  formulations  for  problems  of  large  plastic 
flow,  suggest  a mixed  formulation  based  on  the  functional: 


r2  >> 
'kkekk  “ ¥ kk' 


VN 

• TijV,i>u(k,j)  + \ Tijuk,iuk.j  - p»Bi“i}dv  - S Vids  <88> 


where  e!,  = u’  and  o'  is  related  to  e'  through  an  equation  of  the  type 

of  Eq.  (79a).  It  is  worth  noting  that  the  above  formulation,  Eq.  (88),  is 

analogous  to  the  present  formulation  given  in  Eq.  (87),  except  for  the  fact 

that,  whereas  appears  as  a variable  in  Eq.  (87),  appears  in  Eq.  (88). 

It  is  interesting  to  observe  that  the  procedure  based  on  Eq.  (88)  ceases  to 

be  valid  in  the  limit  of  precise  incompressibility.  Moreover,  in  the  discrete 

(finite-element)  version  of  the  functional  corresponding  to  Eq.  (88)  (when 

0 0 

appropriate  discrete  approximations  for  u^  and  e ^ are  introduced), 

NAGTEGAAL,  et  al.,  (1974)  proceed  to  eliminate  as  a variable  at  the 
element  level  and  introduce  a modified  definition  for  the  strain  energy 
density  functional,  W.  The  rigorous  theoretical  validity  of  the  modified 
discrete  functional,  as  a variational  basis  for  obtaining  discretized 
equilibrium  equations,  appears  somewhat  questionable. 

We  note  that  the  above  discussed  difficulties  with  the  incompressibility 
constraint  are  somewhat  easier  to  handle  in  the  case  of  assumed  stress  finite 
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element  methods  based  on  a complementary  rate  principle  of  type  given  in  Eqs. 

(73  and  74).  [For  a treatment  of  incompressibility,  using  assumed  stress 
finite  element  methods,  see  for  instance,  the  works  of  TONG  (1969),  and 
PIAN  and  LEE  (1976)  in  linear  elastic  infenitesual  deformation  cases  and 
that  of  MURAKAWA  and  ATLURI  (1978b)  in  finite  elasticity  problems]. 

3.2:  Rate  Variational  Principles  in  Total  Lagrangean  Formulation: 

In  the  numerical  solution  of  certain  problems  such  as,  for  instance, 

plates  and  shells,  it  may  be  preferrable  to  use  rate  formulations  wherein  all  the 

variables  in  each  subsequent  increment  are  referred  to  a fixed  Lagrangean 

or  Total  Lagrangean  (TL)  frame.  Thus  in  the  TL  formulation,  the  initial 

configuration  C , with  coordinates  x. , is  used  to  refer  all  the  state 
o 1 

variables  in  each  of  the  subsequent  configurations.  Let  s'  and£'  be  the 
rates  of  2nC*  and  lSt  Piola-Kirchhoff  stresses,  in  going  from  to 
which  stress  rates  are  referred  to  and  measured  per  unit  area  in  the  initial 
configuration  Cq.  Let  V°  be  the  gradient  operator  in  the  coordinates  in 
Cq,  and  set  u be  the  rate  of  displacement  from  the  current  state.  Then  the 
Total  Lagrangean  strain-rate,  jS',  is  given  by, 

E'  = \ [V°u  + (V°u)T  -I-  (V°u)  . (V°uN)T  + (V°uN)  . (V °u ) T]  (89) 

N 

where  u is  the  displacement  at  as  measured  from  Cq.  It  is  seen  that 
the  TL  and  UL  strain-rates  are  related  by: 

NT  ’ N 

E'  = (F  ) . e . F . (90) 

/v  Vv  /v  /v  x 

where  + V°uN)^  and  £ is  defined  in  Eq.  (33).  Likewise,  if  £'  [=(V°u_)^] 

is  the  TL  rate  of  displacement  gradient,  it  is  related  to  the  UL  rate  e by 

e'  =e  . / (91) 
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Also,  from  Eqs.  (3  and  2),  respectively,  the  relations  between  the  TL 
rates  jl'  and  t?  and  UL  rates  ^s  and.^t,,  respectively,  can  be  derived  as: 

a'  = JN(FN)_1s  (FN)‘T 


f .jVr'i.  v . fot  + s" . e'T 

<V  A/  /V  <V  A*  /N/ 


(93a, b) 


N N 

where  J is  the  value  of  the  determinant  of  the  matrix  [y  . .]  . Finally, 

1 » J 

the  TL  rate  of  Jaumann  stress,  ^r',  is  seen,  from  Eqs. (4)  and  (5),  to  be 
related  to  t'  and  s'  by: 

A/ 

£'  = 2 ^ • Q1'  + «'T  • t”1  + t’  . aN  + am  . t,T]  (94) 
= \ [s'  . (I  + hN)  + Q + hN)  . s'  + sN  . h'  + V . aN] 

2 ~ y**  <v  ' )V  ' A/  A«  A«  <v  J 


N N 

whereat  and  ^s  are,  respectively,  the  first  and  second  Piola-Kirchhof f 

stress  tensors  in  C.T  as  referred  to  and  measured  per  unit  area  in  C : 

N O 

N M 

h is  the  engineering  strain  tensor  in  C referred  to  C ; and  a and  a' 

**  N o ^ ^ 

N 

are  rotation  tensors,  such  that  (a,  + a')  is  an  orthogonal  tensor,  and 
these  are  found  from  the  application  of  the  polar -decomposition  theorem  as: 

(V°/)T  = aN  . Q + hN)  and  (V°U)T  = *’  . (I  + hN)  + <*N  . h’ 


Now  we  consider  the  question  of  the  forms  of  rate  potentials  for 

^s',  ^t',  and^r'.  First,  we  note  that  if  a rate  potential  for  s of  the  form 

of  Eq.  (17)  exists,  then,  inview  of  Eqs.  (90)  and  (92),  a rate  potential 

say  W',  can  also  be  shown  to  exist  forks'.  Specifically,  let  the  poten- 
• • 

tlal  W for  s Eq.  (17)  be  written  as: 


W = 2 MijkXeijekX;  mA 


(97a, b) 


f 


Where,  the  tensor  can  be  expressed  interns  of  L^^  and  the  relevant 

plasticity  parameters  through  Eqs.  (18  and  16). 

Then,  in  view  of  Eqs.  (90)  and  (92),  the  rate  potential  W'  forks'  can 
be  written  as: 


where 


W'  = 2MIjk*  EIj  EU> 

N 5xi 

MijkX  J Mrstp  . N 


SW'/dE'  = s' 


dx. 

dx, 

k 

5xe 

„ N 

„ N 

„ N 

9ys 

9yt 

9yp 

(98a, b) 
(99) 


Likewise,  inview  of  relation  (93b)  it  can  be  seen  that  a rate  potential  U' 
exists  for^t',  for  the  case  of  an  classical  elasto-plastic  material,  where, 

U'  - W + x sN:(e'T  . e');  SU'/Se'1  = t'  (100a,b) 

r>j  Vs/  rss  /-s/  / s/ 

In  writing  the  right  hand  side  of  Eq.  (100a),  the  relation  jS1  = 1/2 
T N NT 

(e1  . + J?  . ^g')  may  be  used.  Likewise,  inview  of  Eq.  (95),  it  can 

also  be  seen  that  a rate  potential  Q'  for  the  TL  rate  of  Jaumann  stress, 
r'  also  exists,  where, 

Q'  « W +|sN:  (h'  . h');  9Q’/ah'  = r'  (101a, b) 

Again,  in  writing  the  right  hand  side  of  Eq.  (101a),  the  relation,  11'  = 

1/2  [h'  . (I  + hN)  + (I  + hN)  . h']  may  be  used. 

/V  Vs/  /V  /V  IS/  /V 

Now  we  consider  the  rate  form  of  the  field  equations  and  boundary 
conditions.  Considering  the  rates  of  Eqs.  (6-12c)  these  can  be  written  as: 

In  terms  of  s';  E';  and  u': 

/s/  — 


(LMB>-  V°  . |<sN  . *'T  + s'  • (^N)Tj-  + P°B'  = 0 

,.T 


(AMB)- 


s = s 


(102) 

(103) 
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or,  equivalently, 

(AMB)-*  h'  . / . aN  + QjN  + I)  . <t'  . + tN  . a')  (109) 

= symmetric 

(Compatibility)-*  ^e'  = (V°i)T  (110) 

or,  equivalently,  e'  - a'  . (I  + hN)  + . h'  = (V°u)T  (111) 

(TBC)-*  n . t'  = t'  = t'  at  S (112) 

o 

* -- 

(DBC)-  u = u at  Su  (113) 

o 
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3.2.1:  General  Variational  Principles  in  TL  Rate  Form: 
3 . 2 . 1 . 1 : In  terms  of  j ^ ; and  u 


Using  procedures  analogous  to  those  leading  to  Eq.  (42)  of  the  UL 

case,  a general  rate  variational  principle  governing  Eqs.  (102  to  106,  and 

* 

98b)  can  be  stated  through  the  condition  of  stationarity  of  the  functional 


and  a)  - J <£’>  - pV  • “ + i/’Cs'1  • £') 

:U'  - 1/2  (£•  + e'T  + a'T  • / + • a')]}d 


- s 


(u  - u)  ds 


(114) 


o o 

o * T * 

where  /e'  = (V  u)  ; W is  as  defined  in  Eq.  (98a),  and  t/  is  as  defined  in 
Eq.  (105).  Once  again,  a functional  to  check  the  true  satisfaction  of 
the  fully  nonlinear  field  equations  in  can  be  derived,  but  is  omitted 
here  [see  MURAKAWA  (1978), for  instance,  for  further  details].  If  one 
eliminates  E'  and  s'  from  Eq.  (114),  by  a priori  satisfying  Eqs.  (98b,  104, 
and  106),  one  obtains  a potential  energy  rate  principle  with  the  associated 
functional: 


tt2(u)  = J [W’(u)  - p °B ' . u + 1/2  sN:  (e,T.  e')]dV 
P v 

o 

- J I'  • u ds.  (115) 


o 

* T 

where  e'  - (v u)  . Likewise,  by  inverting  Eq.  98b,  one  can  achieve  a contact 


transformation. 


W 


(116) 


*This  functional  can  be  modified,  in  a manner  analogous  to  that  leading  to 
Eqs.  (82  and  85)  respectively,  to  derive  rate  principles  which  can  be  used 
to  treat  cases  of  near  or  precise  incompressibility. 


11 


such  that 


e'  = asVag1 


(117) 


Using  Eq.  (116),  one  may  eliminate  E'  as  a variable  from  Eq.  (114)  and  obtain 

2 

a functional,  say  tt^^s')  corresponding  to  a HELL INGER -RE I S SNER  type 
rate  principle. 

In  general,  a complementary  energy  rate  principle  may  be  derived  from 
Eq.  (114)  by  eliminating  E 1 from  Eq.  (114)  using  Eq.  (116);  and  by  satisfying 
both  the  equations  of  LMB  and  AMB,  Eqs.  (102)  and  (103)  respectively,  a 
priori.  When  this  is  done,  one  can  formally  obtain  a complementary  energy 
rate  functional, 

^(Hvg')  =J  {s’*(s')  + i,sN:[(V°u)  . (V°u)T]}dV 

V J 

o 


u ds 


(118) 


Thus,  as  in  the  UL  rate  formulation,  even  in  the  TL  formulation,  both  s* 

J 

and  ii  appear  as  variables  in  the  complementary  principle;  however  there  is 

a significant  difference  between  the  two  cases  from  the  point  of  view  of 

application.  In  the  TL  rate  formulation,  the  AMB  condition,  Eq.  (103),  is 

quite  simple  to  be  satisfied,  provided  the  chosen  is  symmetric.  However, 

in  LMB  condition,  Eq.  (102),  both  the  stress  raters',  and  the  displacement 

gradient  rate^e'  are  involved;  moreover,  there  is  a strong  coupling  between 

N 

s'  and  the  currently  known  functions,  F (x) . Thus  the  admissible  stress 
field^s',  to  be  used  in  a complementary  energy  rate  principle,  if  one  were 
contemplated  based  on  £ ' , must  represent  a solution  to  the  set  of  partial 
differential  equations,  Eq.  (102),  with  variable  coefficients.  While  this 
may  mathematically  not  be  impossible,  it  defeats  the  very  purpose  of  a 
variational  principle  forming  the  basis  of  a simple  numerical  method  such 


■ 
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as  a finite  element  method.  Thus  the  rate  complementary  energy  principle 
based  on  j>',  for  finite  strain  plasticity  analysis,  does  not  appear  to  be 
practically  useful. 


3. 2. 1.2:  Based  on  ~t,' ; e';  and  u: 

n*  o* ~ 

Likewise,  a general  variational  principle  governing  Eqs.  (107,  108,  110, 
112,  113  and  100b)  can  be  shown  to  be  governed  by  the  stationarity  condition 
of  the  functional: 


~ S>  = I {U'^'>  ‘ Po5-’  * - 


^ t'  • i ds  > | _t'  . (u1  - u ' ) ds 

a u 


(119) 


o o 

where  U'  is  defined  through  Eq.  (100a);  and  Jt'  is  defined  in  Eq.  (112).  As 
in  the  UL  rate  case,  because  of  the  special  structure  of  U'  as  given  in  Eq. 
(100a),  it  can  be  shown  that^t',  as  derived  from  U'  through  Eq.  (100b), 
identically  satisfies  the  AMB  condition,  Eq.  (108). 

If  from  Eq.  (119)  one  eliminates^'  and  ^t'  as  variables,  by  a priori  ' 
satisfying  Eqs.  (100b,  110,  and  113),  one  can  derive  a rate  form  of  a 
potential  energy  functional : 


TTp(u)  = J |U'(u)  - p°B'  • u]"  dv  - J Jt * . u ds  (120) 

Vo  So 

o 

the  stationarity  of  which  functional  leads  to  Eqs.  (107,  108,  and  112)  as 
its  Euler  equations  and  natural  b.c.  However,  the  above  variational  prin- 
ciple can  be  seen  to  identical  to  that  in  Eq.  (115),  because  of  Eq.  (100a) . 

By  inverting  Eq.  (100b),  one  may,  under  certain  conditions  analogous 
to  those  discussed  in  the  UL  rate  case,  achieve  a contact  transformation. 
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U'^e')  - t'T:e'  = -T'*(V) 


Such  that  dT'*/5t'  = e,T 


(121) 

(122) 


However,  analogous  to  the  situation  discussed  earlier  in  connection  with  the 

UL  rate  formulation,  the  AMB  condition,  Eq.  (108),  cannot  be  verified  to 

be  embedded  in  the  structure  of  T'  (t')  as  obtained  from  Eq.  (121).  Thus 

the  HELLINGER-REISSNER  type  principle  in  terms  of  ^t1  and  u [derivable  by 

using  Eq.  (121)  in  Eq.  (119)],  or  the  complementary  energy  rate  principle 

in  terms  of  <t'  alone  [derivable  by  satisfying  conditions  of  Eqs.  (121,  107, 

and  112)  a priori  in  Eq.  (119)]  cease  to  be  rational  principles;  since,  the 

* 

AMB  condition  for  t'  is  neither  embedded  in  the  structure  of  T'  nor  does 

A/ 

it  follow  as  an  Euler  equation  from  these  principles. 

3 .2. 1 .3:  In  terms  of  a ' ; £ 1 ; u ' ; and  ,r 1 (J: 1 ; g 1 ) 

Once  again,  to  avoid  the  above  difficulties  in  formulating  consistent 
complementary  energy  and  Hellinger-Reissner  type  rate  variational  principles, 
we  transform  the  general  variational  principle  associated  with  Eq.  (119) 
into  one  involving^';  J)';  and  u as  variables.  To  this  end  we  first 
note  from  Eqs.  (100a)  and  (101a), 


U'  = Q'  - 1/2  sN:  (h'  . h')  + 1/2  sN:  (e,T  . e’) 

^ /V  >V  /W  A/  V*  /V 


Upon  making  use  of  the  relations. 


(123) 


rw  /v 


a + hN)  + oN  . h' 

)v«  A# 


(124) 


and  the  orthogonality  condition  and  its  rate  form; 


NT  N T ,T  . N NT 
a .cr=I;  a . cr  = -or 

r*  a^  a/  a# 


Z 


(125) 
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N N NT 

and  the  relation^  = js  . F , one  can,  through  relatively  straightforward 
algebra,  reduce  Eq.  (123)  to 

U’  = Q’  - tNT:  ( a ’ . h’)  - 1/2  tOT:  [o'  . c*NT  . ex'  . (I  + hN)]  (126) 

^ f>*  r*  e*t  f+J  (V  rw 

Upon  using  Eq.  (126)  and  (124),  Eq.  (119)  can  be  rewritten  as: 

«'s  £')  {q‘ <>’)  - pV  • « 

o 

+ t,T:[(V°u)T  - g'  V (I  + tiN)  - aN  . h ’] 

/V  — — >w  sw  /V  J 

-tNT:(Dr'  . h')  - 1/2  t^O'  . . a'  . a + hN)]}dV 

/v  /V  <v  /v  rw  rv  >w  /w  J 

- r t'  . u ds  - r t'  . (u  - u)  ds  (127) 

V_  **  CL  ““ 


Where  Q'  (V)  is  the  potential  for as  defined  in  Eq.  (101a)  and ^r'  is  to 
be  related  to ^t'  and  a'  through  Eq.  (94).  Noting  that  the  variations  6o' 
are  required  to  satisfy  the  property  of  Eq.  (125b),  viz.,  that  . at' 
is  skew  symmetric,  it  can  be  easily  shown  that  the  stationarity  condition  of 
the  above  functional  yields  the  Euler  equations  and  natural  b.c:  (i)  the  LMB 
condition,  Eq.  (107);  (ii)  the  AMB  condition,  Eq.  (109);  (iii)  the  compatibility 
condition,  Eq.  (Ill);  (iv)  the  rate  constitutive  law,  Eq.  (10lb);  (v)  the  TBC, 
Eq.  (112),  and  (vi)  the  DBC  Eq.  (113). 

By  inverting  (even  if  numerically)  the  relation  of  Eq.  (101b),  one 
can  obtain  a contact  transformation. 


Q'  - r':h'  = -R'  (r1) 


Q'  - 1/2  [t'  . aN  + fcN  • a'  + oNT  . t,T  + o'T  . tNT]:h'  = -R  '*<£') 

***  ~ rw  M iv  /w  ^ 


(128) 
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such  that 


(129) 


*R'*/dr'  - h' 

By  using  Eq.  (128),  one  can  eliminate  h'  as  a variable  from  Eq.  (127) 

2 

and  obtain  a functional  rr,^(u';  g';  't')  corresponding  to  a HELLINGER-REISSNER 
type  variational  principle. 

Finally,  by  requiring  t'  to  satisfy  only  the  linear  LMB  condition,  Eq. 

(107),  and  the  TBC,  Eq.  (112),  one  can  eliminate  u1  as  a variable  from 

“ 

2 

TT^  and  obtain  a TL  rate  complementary  energy  functional, 

"o  <£'%'>  - J {-R'V>  • a + iN>] 

Vo 


.NT  r , 


NT 

or 


a 


(I  + hN)]}dV 


+ [ t ' . u ds 
u 

o 


(130) 


Wherein  it  is  implied  thatjr'  is  related  to  t'  and  through  Eq.  (94). 

Noting  that  the  variations  6^t'  are  now  subject  to  the  constraints  that 
V°  . 6t'  = 0 in  V and  n . 6t*  = 0 at  S ; and  that  the  variation  6a'  are 

“ yw.  O “ ^ O 

o 

NT 

subject  to  the  constraint  that  a . 6g ' is  skew- symmetric,  it  can  be  shown 
easily  that  the  condition  of  vanishing  of  the  first  variation  of  the  above 

_ * m »r 

functional  leads  to:  (i)  compatibility  condition,  (VM*)1  = or'  . (1+  hN)  + 
&N  • ij';  (u)  the  A*®  condition,  Eq.  (109);  (iii)  the  DBC,  Eq.  (113). 

Once  again,  in  as  much  as  the  AMB  condition  follows  unambiguously  as 
in  Euler  equation;  and  the  admissible^'  is  required  to  satisfy  only  the 
uncoupled  linear  LMB  condition,  Eq.  (107)  [which  can  be  satisfied  easily 
through  first-order  stress  functions,  as  ^t'  = V°  x JF  + ^t'p;  arid  jt'P  is  any 
particular  solution  such  that  V°  . (t'p=  -p°B'],  and  the  TBC,  Eq.  (12),  the 


< 
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TL  rate  complementary  energy  principle  of  Eq.  (130),  is  the  most  consistent 
and  useful  principle  for  purposes  of  engineering  application. 

As  in  Eq.  (75)  for  the  UL  rate  case,  in  the  application  of  Eq.  (130) 
also  to  a finite  element  assemblage,  the  interelement  traction  reciprocity 
condition,  (n  . t')+  + (n  . t')  at  p can  be  relaxed  a priori,  and  in- 
troduced  as  a constraint  condition  in  a modified  complementary  energy 
principle  through  Lagrange  Multipliers  ^ at  p^.  This  modified  principle 
is  stated  through  the  functional: 

tt^  ( t ' ; a' ; u = > f -f-R'  (r')  - t,T:[a'  . (I  + hN)] 
mc  ~ ~ — p L. i J l ~ ~ ~ ~ ~ J 

m V 


1 .NT  r | NT  I . T i 1N.-1  I Jtt 

'2^  :[£  £ • Z ■ +iJ  >^/dv 

T r t-  . u ds  + V j (a  . t-)  .z  dp 

m S m p 


(131) 


As  in  the  UL  rate  case  [see  discussion  following  Eq.  (75)]  the 
functional  in  Eq.  (31)  can  be  used  to  develop  a finite  element  "stiffness 
matrix"  method  and  Eq.  (130)  can  be  used  to  develop  finite  element  "flexibility 
matrix"  approach.  The  author's  work  in  progress  in  this  regard,  along  with 
analysis  of  certain  metal  forming  problems  is  the  subject  of  a forthcoming 


paper. 

Finally  it  should  be  remarked  that  eventhough  the  development  of  Eq.  (130), 
as  a basis  of  a TL  rate  complementary  energy  principle  for  finite  strain 
elasto-plasticity  was  based  on  independent  considerations,  the  result 
is  analogous  to  the  principle  derived  by  FRAE1JS  DE  VEUBEKE  (1972).  However, 
FREIJS  DE  VEUBEKE 'S  (1972)  principle  governs  the  total  (as  opposed  to  rates) 
deformations  of  a compressible  elastic  solid.  It  is  also  noted  that  a TL 
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rate  principle  equivalent  to  that  of  FRAEIJS  DE  VEUBEKE  (1972)  was  develooed 
and  used,  in  the  context  of  an  assumed  stress-rate  finite  element  method, 
to  solve  certain  finite  strain  problems  for  nonlinear  elastic,  compressible 
as  well  as  incompressible  solids  by  MURAKAWA  and  ATLURI  (1978a, b). 

Finally  we  remark  that  the  TL  rate  complementary  energy  principle  for 
elastic-plastic  solids  given  presently  in  Eq.  (130)  differs  slightly  from 
that  for  nonlinear  elastic  solids  given  eariler  by  MURAKAWA  and  ATLURI 
(1978a),  this  difference  is  in  the  third  term  in  the  volume  integrand  on 
the  right  hand  side  of  Eq.  (130)j  instead  of  the  term  appearing  in  Eq.  (130), 
the  term  - ^t^ : [<* ' . (J[  + h1*)]  appears  in  the  paper  by  MURAKAWA  and  ATLURI 
(1978a).  The  effect  of  this  is;  whereas  the  exact  rate  condition  of  AMB, 

Eq.  (109)  becomes  an  Euler  equation  of  the  principle  in  Eq.  (130),  the  AMB 
condition  which  follows  as  an  Euler  Equation  for  the  principle  given  in 
MURAKAWA  and  ATLURI  (1978a)  is  that  (I  + hN)'  . t'  . c*N  + h'  . tN  . <*N  is 
symmetric  [which  condition  represents  only  an  approximation  to  the  exact 
condition,  Eq.  (109)].  However,  the  iterative  correction  procedures  to 
check  the  satisfaction  of  the  fully  nonlinear  AMB  condition  of  Eq.  (11a)  at 
the  end  of  each  increment  were  employed  by  MURAKAWA  and  ATLURI  (1978a)  to 
correct  the  above  approximation  in  the  rate  AMB  condition.  In  this  sense, 

the  principle  currently  stated  through  Eq.  (130)  [which  is  equally  applicable 

* 

to  nonlinear  elastic  solids,  when  the  potential  R'  Q:  ) is  appropriately 
defined]  is  the  most  consistent  TL  rate  complementary  energy  principle  for 
finite  strain  analysis  of  elastic  as  well  as  elastic-plastic  solids. 
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